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Abstract. In order to estimate the masses of the com- 
pact, young star clusters in the merging galaxy pair, NGC 
4038/4039 ("the Antennae"), we have obtained medium 
and high resolution spectroscopy using ISAAC on VLT- 
UTl and UVES on VLT-UT2 of five such clusters. The 
velocity dispersions were estimated using the stellar ab- 
sorption features of CO at 2.29 /im and metal absorption 
lines at around 8500 A, including lines of the Calcium 
Triplet. The size scales and light profiles were measured 
from HST images. From these data and assuming Virial 
equilibrium, we estimated the masses of five clusters. The 
resulting masses range from 6.5 x 10^ to 4.7 x 10^ Mq. 
These masses are large, factor of a few to more than 10 
larger than the typical mass of a globular cluster in the 
Milky Way. The mass-to-light ratios for these clusters in 
the V- and K-bands in comparison with stellar synthesis 
models suggest that to first order the IMF slopes are ap- 
proximately consistent with Salpeter for a mass range of 
0.1 to 100 Mq. However, the clusters show a significant 
range of possible IMF slopes or lower mass cut-offs and 
that these variations may correlate with the interstellar 
environment of the cluster. Comparison with the results of 
Fokker-Planck simulations of compact clusters with prop- 
erties similar to the clusters studied here, suggest that 
they are likely to be long-lived and may lose a substantial 
fraction of their total mass. This mass loss would make 
the star clusters obtain masses which are comparable to 
the typical mass of a globular cluster. 



young star clusters which formed during the merging 
process ( Holtzman et al. 1992 



Meurer et al. 1995 



Whitmore et al. 1993 



Whitmore et al. 1997 



Whitmore et al. 1999, Zepf et al. 1999). The overaU 



properties of these clusters suggest that they could be 
the progenitors of the globular cluster populations seen 
in normal nearby elliptical and spiral galaxies (e.g., 
^epf fc Ashman 1993 ). Such an hypothesis also has the 
attractive implication that if ellipticals formed through 
the merger of two large spiral galaxies (as suggested in the 
popular hierarchical merging model of galaxy formation 
e.g., Kauffmann, White, and Guiderdoni 1993), these 
young clusters might evolve into the red (supposedly 
metal rich) part of the globular cluster population of 
elliptical galaxies when the merger is complete (e.g. 
^chweizer 2001 ). However, this hypothesis needs to 
be tested by determining the characteristics of both, 
individual young clusters and the cluster population as a 
whole. 

Globular clusters have a typical mass of 1-2x10^ 
Mq and a mass function which is log-normal 
pJarris 1991 ). The population of clusters in 
4038/4039 ("the Antennae"), however, has a power 
law luminosity function, and the same shape is also 



(,e-g-, 
NGC 



suggested for the mass function ( Whitmore et al. 1999 , 
^hang fc Fall 1999 ). Masses determined from photometric 
data are as large as a few x 10^ Mq for some of the clus 



ters (Zhang fc FaU 1999, Mengel et al. 2001) and the de 
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termined ages span a large range (Whitmore et al. 1999 



dynamical masses - NGC Mengel et al. 2001) 



1. Introduction 

During the last years, many interacting and merging 
galaxies were discovered to host large numbers of 
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At first glance, the population of globular clusters in 
the Milky Way and those in the Antennae seem to have 
vastly different ensemble characteristics. For example, the 
most massive clusters in the Antennae are at least a fac- 
tor of few more massive than those in the Milky Way. 
However, given the large number of clusters formed in the 
Antennae, it seems reasonable that the mass function is 
sampled up to a high upper mass, and moreover mass loss 
during the evolution of the young clusters can be expected 
to play a role. Evolution over a Hubble time might convert 
the power law cluster mass function into a log-normal one, 
if lower mass clusters were dissolved preferentially during 
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the evolution. Dynamical models like those, for example, 
of Chernoff & Weinberg (1990) and Takahashi & Portegies 
Zwart (2000) provide theoretical support for the necessary 
differential evolution from a power-law mass distribution 
function into a log-normal distribution. 

Dynamical cluster masses are less model dependent 
than those determined from photometry only, or, at the 
very least, have a different set of dependencies. The dy- 
namical mass is derived from the stellar velocity disper- 
sion in combination with the cluster size and light pro- 
file. While the photometric mass estimates depend on the 
model assumed for the star formation parameters (time- 
scale, age, IMF slope, lower and upper limiting masses, 
metallicity, etc), the dynamical mass estimates rely only 
on the validity of the Virial equilibrium (i.e., the potential 
is due to the collective gravitational effects of individual 
stars - self-gravitating - and is changing only slowly with 
time), and the assumed constancy of the M/L ratio within 
the cluster. The comparison of photometric and dynam- 
ical cluster masses (cluster M/L) constrains the slope of 
the mass function and may reveal the presence or absence 
of low mass stars. The fraction of low mass stars influences 
the survival probability of the cluster during a few Gyrs of 
evolution, as clusters rich in low-mass stars are less prone 
to destruction (e.g., Takahashi & Portegies Zwart 2000). 

More specifically, we have used high spectral resolu- 
tion spectroscopy conducted at the ESO VLT to esti- 
mate the stellar radial velocity dispersion tr, and high 
spatial resolution imaging from archival HST images 



(Whitmore et al. 1999) to estimate the size scales (e.g. 



the projected half-light radius r/jp). This results in Mdyn 
for individual clusters which can be estimated using the 
equation: 

M = ^ (1) 

Where 77 is a constant that depends on the distribution 
of the stellar density with radius, the mass-to-light ratio 
as a function of radius, etc. The only assumption that 
goes into Virial relations of this form is that the cluster is 
self-gravitating and is roughly in equilibrium over many 
crossing times (i.e., it is not rapidly collapsing or expand- 
ing). The calculations of Spitzer (1987) indicate that for 
a wide range of models, 77 w 10. Under the assumption of 
isotropic orbits, i] differs from 3 (the pure Virial assump- 
tion) since the "gravitational radius" must be scaled to 
the projected half-light radius (or some measurable scale) 
and because the measurement of the velocity dispersion is 
not the central velocity dispersion but is a measurement 
which is weighted by the light profile. So differing mass 
profiles have different values of 77 appropriate for estimat- 
ing their dynamical masses. For example, globular clusters 
in the Milky Way exhibit a range of concentration param- 
eters (logarithm of the ratio of the tidal radius to the core 
radius) in a King model (King 1966) of approximately 0.5 
to 2.5 (e.g., Harris 1991). Over this range of concentra- 
tions, rj varies from about 9.7 to 5.6. Hence, accurately 



measuring the light profiles (and of course the assumption 
that the light profile traces the mass profile) is crucial in 
estimating the masses. 

Recent attempts at estimating the dynamical masses 
of similar clusters in other galaxies has met with 
Notably, ' 



some success 



Ho 



Ho & Filippenko 1996a and 
Filippenko 1996b| were able to estimate the masses 



of two luminous young star clusters in the blue compact 
dwarf starburst galaxy NGC 1705 (1705-1) and NGC 1569 
(1569A). They measured velocity dispersions for 1705-1 
and 1569A of 11.4 km s~^ and 15.7 km s~^ respectively. 
Under the assumption that rj = 10, Sternberg (1998) es- 
timated that these clusters have masses of about 2.7 x 
10^ Mq for 1705-1 and 1.1 x 10^ Mf^ for 1569A. More 







recently, Smith & Gallagher (2001) have estimated the dy- 
namical mass of the luminous star-cluster in M82, M82-F. 
Assuming 77=10, they find a mass of 1.2±0.1x 10^ Mq. 
This mass estimate implies that the mass-to-light ratio of 
M82-F is very extreme and requires either a very flat mass 
function slope or a lack of low mass stars for a Salpeter- 
like mass function slope. Interestingly, this small sample 
of clusters seems to require a range of IMFs to explain 
their mass-to- light ratios. 

The purpose of the present paper is to derive masses 
for a small sample of young compact clusters in the Anten- 
nae galaxies - the nearest merger - and to use these mass 
estimates to constrain the functional form of the IMF in 
these clusters. Since it is important that the velocity dis- 
persion be determined from the stellar component of the 
cluster, we have undertaken a program to observe extin- 
guished clusters in the K-band, to take advantage of the 
strong CO absorption band-heads beyond 2.29 ^m, and 
unextinguished clusters in the optical, to take advantage 
of the strong absorption of the Calcium Triplet (CaT) 
and other metal lines around 8500A. These features are 
strong in atmospheres of red supergiants which would be 
expected in large numbers in clusters with young ages {~6- 
50 Myrs) and thus particularly well-suited for estimating 
the masses of clusters like those in the Antennae. Given 
that previous investigators (e.g.. Ho & Filippenko 1996a, 
Filippenko 1996b| ) have measured velocity disper 



Ho 



sions around 15 km s^^ suggests that to conduct such 
measurements would require a minimum resolution of 
R«9000. Such resolutions are now available using ISAAC 
on VLT-UTl which with its narrowest slit delivers 
R«9000 and UVES on VLT-UT2 which delivers a reso- 
lution of «38,000 with a 1" wide slit. 



2. Observations 

2.1. Moderate Spectral Resolution K-band Spectroscopy 

Spectroscopy in the K-band was performed with ISAAC at 
VLT-UTl on April 17 2000. The spectrograph was config- 
ured to have a 0'/3 slit and the medium resolution grating 
was tilted such that the central wavelength was 2.31/im 
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Fig. 1. ISAAC Ks-band image of NGC 4038/4039 (taken 
during a run in April 2001) with labehng of the clusters 
that were observed using ISAAC and UVES spectroscopy 



with a total wavelength coverage from 2.25 to 2.37/im. 
This was sufficient to include the ^^CO (2-0) and (3-1) ab- 
sorption bands and at a resolution (A/ A A) of about 9000. 
Observations of late-type supergiant stars were taken so 
that they could be used as templates for the determina- 
tion of the velocity dispersions^. Various templates were 
observed, however, due to incorrect classifications in the 
literature, only one of them actually was an M-type super- 
giant. More template stars were observed by Linda Tac- 
coni during an observing run on August 15 2000 using the 
same instrumental set-up. Observations were performed 
by nodding along the slit and dithering the source posi- 
tion from one exposure to the next. Atmospheric calibra- 
tor stars were observed several times during the night. 

The target clusters were selected using the results 
from NTT-SOFI imaging observations performed in May 
1999 (Mengel, Ph. D. Thesis, der Ludwig-Maximilians- 
Universitat Miinchen, 2001; Mengel et al. 2001, in prepara- 
tion). They included narrow band imaging in the CO (2-0) 
band-head and revealed those clusters which are at an age 
where the near-infrared emission is dominated by red su- 
pergiants. The clusters with the highest equivalent width 
in the CO band-heads, highest luminosity, and high ex- 
tinction were selected as the ISAAC targets (for the lo- 
cation of all observed clusters see Fig. |^). Two target 
clusters were positioned in the first slit position. These 



were two clusters with a separation of ^ 5" located in the 
region where the galaxy discs apparently overlap ("the 
overlap region" ) . One of them has an optical counter- 
part, [W99]15, the other one is detectable only in the I- 
band, therefore denoted as a "Very Red Object (VRO)" 
by Whitmorc & Schwcizer 1995 (with the designation: 
[WS95]355). Their detectability in at least one band from 
the HST imaging program of Whitmorc et al. 1999| was a 



necessary additional constraint, since the compactness of 
each cluster means that their radii must be measured from 



^ This research has made use of the SIMBAD database, op- 
erated at CDS, Strasbourg, France 



HST data (see §3.4). The second slit position included only 
one bright cluster ([W99]2), which has a relatively low ex- 
tinction and was also observed with UVES. Due to the 
extreme narrowness of the slit, object drifts made several 
re-acquisitions necessary, especially when the object was 
transiting. The ISAAC observations are summarized in 
Table 2. 

2.2. High- Resolution Optical Spectroscopy 

High-resolution optical echelle spectroscopy was per- 
formed using UVES at VLT-UT2 on the night of 18 April 
2000. The instrument was configured with a slit width of 
1", which resulted in a resolution of R«38,000 (depends 
on wavelength). A dichroic was placed in the light path 
allowing the use of both the red and blue arms of the spec- 
trograph. However, in this paper we will discuss only the 
results obtained with data from the red arm. The central 
wavelength of the red arm was shifted to 8400A since part 
of the CaT (at NGC 4038/4039 redshift of z~0.005 located 
roughly at: 8540A, 8585A and 8705A) would have fallen 
right in the small gap between the two CCDs that cover 
the lower and the upper part of the spectra in the red arm. 
The CCDs were read out with a binning of 2x2, 50kHz and 
high gain. For the cluster selection we applied the same 
criteria as for the clusters observed with ISAAC, how- 
ever, for the UVES observations, we selected only those 
with relatively low optical extinction and bright I-band 
magnitudes. To boost the efficiency of the observing, we 
additionally required the clusters to have a nearby clus- 
ter within the slit of UVES (which is only 10" due to the 
echelle design). The clusters observed were [W99]l and 
[W99]43 in one slit position, [W99]10 and [W99]16 in a 
second slit position. Unfortunately, [W99]43 and [W99]10 
turned out to be too faint to extract a spectrum of suffi- 
cient S/N for further analysis. A third slit position covered 
[W99]2, common to ISAAC and UVES. Obtaining two in- 
dependent observations of [W99]2 provides an indepen- 
dent estimate of the velocity dispersion to test the uncer- 
tainties of our measurements and to gauge whether there 
might be some systematic differences in velocity disper- 
sions determined from the CaT in the optical with UVES 
versus that obtained from the CO band-head with ISAAC 
in the near-IR. A range of template supergiants was ob- 
served (early K through late M-type supergiant stars) , but 
also several hot main sequence stars (late O through mid 
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Table 1. Summary of ISAAC spectroscopy. The first slit position included two clusters, the second one cluster. The 
supergiants were observed to provide template spectra for the determination of velocity dispersions. The slit width 
was O'.'S, resulting in a resolution R = A/AA of 9000. Seeing variations are expected to have no effect on the spectral 
resolution given the narrow slit width. The total integration time was split up into several integrations of duration 
Tint- The K-band magnitudes were determined from an ISAAC Ks-band image (aperture diameter 1'.'74). The seeing 
value refers to an average for the monitored V-band seeing during the integration. The stars were observed during two 
observing runs (2000 and 2001), and those marked with LT were kindly observed by Linda Tacconi during another 
run, because the 2000 run had resulted in only one useful supergiant spectrum. 



Position 


Object 


mv 


niK 






Total ON time 


Seeing 


Comments 






[mag] 


[mag] 


[mag] 


M 


[s] 






Slit 1 


[W99]15 


19.4 


15.9 


15.8 


600 


15600 


0'.'7 


Difficulties keeping both 




[WS95]355 


^22.1 


15.7 


15.4 


600 


15600 




targets in slit during 2 integrations 


Slit 2 


[W99]2 


17.0 


14.0 


14.0 


600 


2400s 


o'.'g 





and the following stars: K5Ib (HD183589), K5Ib (HD53177, LT), MOI (CD-44 11324, LT), 

MO-lIab (HD 42475), Mil (BD+0 4030), Mllab (HD 98817), Mllb (HD 163755), M2Iab (SAO 15241), 

M3Ia-Iab (HD103052), M3.5I (BD+13 1212), M5Iab (HP 142154) 



B-type main sequence stars) . Massive main sequence stars 
are expected to dominate the blue spectrum and to con- 
tribute significantly to the flux in the I-band (Bruzual & 
Chariot 1993). Due to the high dispersion and faint sky 
background, integration times could be very long and were 
only limited by the desire to keep the number of cosmic ray 
hits down to a reasonable limit. The UVES observations 
are summarized in Table 2. 

3. Reduction and Analysis 
3.1. Reduction of the near-IR data 

The reduction of ISAAC data was performed in a standard 
way using the IRAF data reduction package^. The reduc- 
tion included dark subtraction and flat-fielding (using a 
normalized flat-field created from dome flat observations) 
on the two-dimensional array. Sky subtraction was per- 
formed pairwise, followed by a rejection of cosmic ray hits 
and bad pixels. The spectra were then corrected for tilt 
and slit curvature by tracing the peak of the stellar spatial 
profile along the dispersion direction and fitting a polyno- 
mial to the function of displacement versus wavelength. 
The corresponding procedure was also applied to the or- 
thogonal direction. Single integrations were combined by 
shifting-and-adding, including a rejection of highest and 
lowest pixels. The rejection could not be applied to ob- 
servations of [WS99]2, because they consisted of only a 
limited number of single frames. The object spectra were 
then extracted from user defined apertures. A linear fit to 
the background (below ?a7% of the peak intensity for all 
clusters) on both sides of the object spectrum was sub- 
tracted. 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



The spectra were wavelength calibrated using obser- 
vations of arc discharge lamps. An atmospheric calibrator 
(B5V) was observed and reduced in the same way as the 
target and used to divide out the atmospheric absorption 
features from the spectra. 

3.2. Reduction of the optical echelle data 

Reduction of the UVES data was performed within the 
echelle spectra reduction environment of IRAF. Reduc- 
tion included bias subtraction, flat-fielding and extracting 
the spectrum from user supplied apertures. A background 
is fit to neighboring regions and subtracted. Also here, the 
background was below »7% of the peak intensity for all 
clusters. This step includes a bad pixel rejection. The full 
moon was fairly close to the Antennae during the integra- 
tions, and its scattered light introduced a solar spectrum 
into the data. However, this contribution is very effectively 
removed during background subtraction. 

Wavelength calibration required the identification of 
many lines in the ThAr-spectra for each of the echelle 
orders which were identified using a line table available 
from the ESO-UVES web pages. A dispersion function 
was then fit and applied to the object data. This data set 
contains each order of the spectrum in a separate channel. 
These separate channels are combined into one final spec- 
trum covering the total wavelength range. This involves 
re-gridding of the wavelength axis into equidistant bins 
and averaging in the overlapping edges of the orders. The 
spectral resolution obtained is R«38,000 across all orders. 

3.3. Estimating the Velocity Dispersions 

For each spectrum (both the optical and near-IR spectra) , 
we estimated the a of the broadening function (assumed to 
be a Gaussian) which best fit the cluster spectrum in the 
following way. The stellar spectrum (the template spec- 
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Table 2. Summary of UVES spectroscopy of science targets (star clusters and template stars). The spectral resolution 
was R~38,000, and differences in the seeing for different objects had no measurable effect on the shape of the absorption 
features in the stellar spectra. 



Position 


Object 


mv 


niK 


mA'o 




Ttot 


Seeing 


Comments 






[mag] 


[mag] 


[mag] 


[s] 


[s] 






Slit 1 


[W99]l 


17.6 


14.7 


14.7 


4800s 


9600s 


0'.'7 


required pausing of integration 




[WS95]43 


20.0 






4800s 


9600s 




due to A.O. software crash 


Slit 2 


[W99]10 


19.0 


15.6 


15.6 


4800s+1800s 


11400s 


0'.'85 


required pausing of integration 




[W99]16 


19.0 


15.5 


15.5 


4800s+1800s 


11400s 




due to A.O. software crash 


Slit 3 


[W9912 


17.0 


14.0 


14.0 


4600s 


4600s 


1'.'2 





And the following stars: M3Iab (HD 303250), M5I (HD 142154), Mllb (HR 6693), Mla-Ib (HR 4064) 
K7IIa (HD 181475), K5Ib (HR 7412), K5II (HR 7873), K3Ib-IICN (HR 6959), K3H (HR 6842) 
K2.5Hb (HR 7604), K2H (HR 6498), K2Ib (BM Sco), K0.5Hb (HR 6392), KOIa (HR 6392) 
B5V (HR 5914), B5V (HR 5839), B2V (HR 6028), B2V (HR 6946) with integration times 
between Is and 20s, and some atmospheric and flux calibrators 



trum) was broadened with Gaussian functions of variable 
tr in velocity space. The resulting set of spectra were then 
compared with the cluster spectrum. The best fit was de- 
termined by evaluating ^ and then search for the min- 
imum of the function x^(vr,cr) using a simplex downhill 
algorithm for the tour through parameter space. 

3.3.1. Unique Character of the near-IR Estimates 

Obviously, the template spectrum has to be a good over- 
all match to the cluster spectrum, otherwise an erroneous 
velocity width will be derived for the cluster. For a star 
cluster that formed '^10 Myrs ago, late K through early 
M supergiants are expected to provide the largest con- 
tribution to the flux at 2.3/im. However, according to 



population syn thesis models (e.g., Lcithcrer et al. 1999 , 
Sternberg 1998), there will be a non- negligible contribu- 
tion to the flux from hot main sequence stars. Since the 
stars are hot (O and B-type stars), this "diluting contin- 
uum" will be an essentially featureless continuum which 
solely decreases the equivalent width of the CO band- 
heads. This has the effect of shifting the apparent domi- 
nating stellar type towards higher effective temperatures. 
Starting out with a template spectrum with weak CO fea- 
tures leads to very low velocity dispersions, the opposite is 
the case if an MSI star (strong band-heads) is used, with 
vast differences in the results (a few km s^^ to up to about 
30 km s^i). 

To first order, the velocity dispersion determined from 
different stellar templates agrees if the CO equivalent 
widths match. This is achieved by adding a (positive or 
negative) continuum to the stellar spectrum, which dilutes 
or enhances the CO features. 

To second order, there are two different indicators that 
can be used to constrain the matching stellar template. 
Both relate to the shape of the CO band-head, due to vari- 
ations in that shape between the different supergiant spec- 
tral classes. The various rotational transitions are resolved 



in the 12CO(2-0) and the i2C0(3-l) band-head. We per- 
formed tests of our analysis technique by creating a sim- 
ulated cluster spectrum with an input velocity dispersion 
(Tin, added noise (usually S/N = 15) and re-determined 
the velocity dispersion (Jout using different template star 
spectra. With the correct template (out — in) , we obtained 
in 30 test runs an average Oout — 14.7±0.7 km s~^ (with 
Gin — 15.0 km s^^). Performing the same test on an input 
spectrum that was diluted by a flat continuum (10%), but 
still fltting with the undiluted template, the velocity dis- 
persion increased (17.9±1.7 km s~^), as expected. Fitting 
other stellar types (MOI, Mil) yields similar results. But 
the mismatch can be diagnosed by: 

— determining a over several different wavelength ranges. 
With the matching template, Oout depends only weakly 
(maximum differences ^ 6%) on the selected wave- 
length range, no matter if only the band-head, only the 
overtones, or both are covered. A template mismatch 
causes differences in the Gout for different wavelength 
ranges of 25% typically, but sometimes more. 

— inspecting the shape of the flt around the tip and the 
rising edge of the band-head, which provides a good 
match only for the correct template 

3.3.2. Unique Character of the Optical Estimates 

The determination of the velocity dispersion for the opti- 
cal echelle data used the same procedure as for the ISAAC 
spectra, relying on the Calcium Triplet around 8500A, but 
also using the Mg absorption feature at 88OOA and other 
weaker metal absorption lines between 8400 and 9000A. 

Examining the spectra in detail, it was obvious that 
they are composed of a mixture of light from supergiants 
and hot main sequence stars. This means that the spectra 
cannot be simply fit by broadening the stellar template, 
but that either the hot star contribution needs to be added 
to the template, as was done in the case of the ISAAC 
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Fig. 2. ISAAC Spectra of the three clusters [WS95]355, [W99]15 and [W99]2 (top to bottom, black) with best fits 
(grey). The bottom plot shows the residuals of the fit to spectrum [W99]15 for the best fit and for two fits with lower 
and higher velocity dispersions. This gives a visual impression of how much the velocity dispersion affects the spectral 
features. 



spectra, or that the contribution needs to be subtracted 
from the cluster spectrum. 

Here, the procedure was different from that applied to 
ISAAC data, because the hot main sequence stars were 
not completely featureless over the analyzed wavelength 
range, but had some Paschen absorption features. Pal6, 
Pais and Pal3 lie slightly red-ward of the CaT wave- 
lengths 8498A, 8542A, and 8662A, respectively. Since the 



widths of these lines are mainly due to line broadening in 
the atmosphere and rotation of the star itself and only an 
insignificant fraction of the width is due to the velocity dis- 
persion of the star cluster, they cannot be used to estimate 
the velocity dispersion. Therefore this weak contribution 
was subtracted from each cluster before the velocity dis- 
persion was estimated. The amount of subtraction was 
estimated by the strengths of the Paschen line absorption 
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Rest wavelength [lam] 

Fig. 3. Displayed are parts of the normalized spectrum 
of cluster [W99]2 (black), together with the fit (red). The 
four pieces were fit separately, and the average value of 
the velocity dispersion is cr = 14.2 km s~^. It is worth 
noting that this is the cluster which was also observed 
with ISAAC, and to compare the results (see Figure ||). 



and the appropriate template for the subtraction was cho- 
sen for its ability to remove the contribution accurately. 
Given the S /N of the spectra and the weakness of the fea- 
tures removed, the final characteristics of the subtracted 
spectrum were not very sensitive to the exact template 
used. 

The cluster [W99]2 merits some specific discussion 
since it has both ISAAC and UVES observations and also 
some peculiarities were noted compared to the other clus- 
ters observed with UVES. Fig. ^ shows fits of an MSI spec- 
trum to pieces of the spectrum of the cluster [W99]2. The 
latter had a 45% contribution of B2V stars subtracted. 
This contribution was determined by interactively sub- 
tracting the spectrum of the B2V star until the Paschen 
absorption features were reduced down to the noise level 
of the cluster spectrum. In addition, we fit the continuum 
of the cluster spectrum using a combination of a B2V and 
an MSI star spectrum, in which the relative contributions 
of the two template stars was the fit parameter. The best 
fitting combination again included a 45% contribution of 
B2V stars. Moreover, the relative contributions of MSI and 
B2V stars in I and K band agree with the determined con- 
tributions to the ISAAC and UVES spectra, respectively, 
for this cluster. 

3.4- Estimating the Sizes of the Observed Clusters 

In order to make an estimate of the mass, it is necessary 
to have a measure of the cluster radius and profile shape. 
Whitmore et al. (1999) found that many of the young star 
clusters in images taken with the Wide Field Planetary 
Camera on HST (WFPC-2) of the Antennae were slightly- 



to- well-resolved and estimated their sizes. Their main aim 
was to determine the average radii of the clusters and to 
compare them to those observed in other galaxies, both 
to the young clusters in other mergers and to globular 
clusters. They found that the mean effective radius of the 
clusters in the Antennae is leff ~ 4 pc. However, Whit- 
more et al. (1999) did not present their results on indi- 
vidual clusters. Brad Whitmore kindly provided us with 
their estimates of the effective radii for the clusters we ob- 
served with ISAAC and UVES. But given the importance 
of size and light profile to the mass determination, we ad- 
ditionally measured them ourselves from the HST I-band 
images. 

Similar to what Whitmore et al. (1999) did to make 
such estimates we reduced archival HST I-band images 
(see Whitmore et al. 1999 for the details concerning the 
images) using the drizzle technique (see the Drizzle Hand- 
book available from STScI). Given the pixel sizes of the 
HST chips (O'.'lOl pix~^ for the Wide Field Camera arrays 
and 0'.'045 pix~^ for the Planetary Camera array), and the 
distance to NGC 40S8/40S9, which implies a scale of 9S pc 
arcsec"^, the factor of two improvement in the sampling of 
the HST point-spread-function (PSF) due to drizzling the 
images is critical in getting a robust estimate of the sizes. 
At the typical cluster size, Tef / ~ 4 pc, the subtended scale 
is only 0'.'043 or about one pixel of the Planetary Camera! 
For this analysis we used the ishape routine described by 



Larsen 1999 implemented in the BAOLAB data reduction 



package developed by Larsen. The ishape routines were 
especially designed to work more efficiently than typical 
deconvolution routines at determining the sizes of slightly 
extended objects, ishape convolves the user-provided PSF 
with the given cluster profile and determines the minimum 
for a range of sizes using a simplex downhill algorithm 
to find the minimum. The PSF we selected for this con- 
volution is one from the Tiny Tim program (Krist 1995) 
which Whitmore et al. (1999) argue persuasively is ap- 
propriate for these data. We fit each cluster with a range 
of provided profiles which include Gaussian p rofiles, Kin g 
profiles of several concentration parameters ( King 1966 ), 
and Moffat profiles of two different concentrations (Moffat 
1967). The program produces a reduced to measure the 
significance of the fit and various images (artificial cluster 
image, the original input image, an image of the residu- 
als between original and fit, and an image of the weights 
applied to each pixel) to judge the quality of the fit. The 
code also provides a handy weighting procedure that en- 
ables the rejection of "outliers" in the radial cluster profile, 
thereby excluding, for example, neighboring stars and/or 
clusters. 



4. Results 

Our results are summarized in Table S. We find a range 
of cluster velocity dispersions from about 9 km s~^ to 
over 20 km s~^. One concern about any such measure- 
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Table 3. The cluster masses as they were derived from the ISAAC and UVES spectra (in column "Inst." indicated by I 
and U, respectively) in comparison with the supergiant template spectrum. The age was derived from the combination 
of W ar7, Wco, and Wg-aT - For cluster [W99]2, it agrees well with the age estimated from the UV spectrum of 7±1 



Myrs Whitmore et al. 1999. The size is the projected half-light radius ihp estimated as described in § p.4| (for the 
cluster marked with °, we could not obtain a satisfactory fit and used the value provided by B. Whitmore). Wco was 
estimated from the ISAAC spectra between rest-frame wavelengths 2.2924/j,m and 2.2977/im, ([W99]2 also from 3D 



integral field spectroscopy), and WcaT from the UVES spectra, according to the definition in Diaz et al. 1989, The 



stellar velocity dispersion a was determined as described in §3.3, for [W99]2 the average value of ISAAC and UVES 
measurements is given (cr=14.0±0.8 and 14.3 ±0.5 km s~^, respectively). Mu^^ is the Virial mass determined from 
equation |l| with the Icr-uncertainties given in column ctm ■ The light-to-niass ratios were derived from the extinction 
corrected magnitudes and a distance modulus to NGC 4038/4039 of 31.41. 



Cluster 


Inst. 


Wco 


WcaT 


Age 


a 




log Mvir 




log Lv/M 


log Lk/M 




I/U 


[A] 


[A] 


[lOVr] 


[km/s] 


[pc] 


[Mq] 


[%] 


[Lq/MqI 


[Lq/Mq] 


[WS95]355 


I 


16.3±0.2 




8.5±0.3 


21.4±0.7 


4.8±0.5" 


6.67 


12 




1.05 


[W99]15 


I 


17.0±0.2 




8.7±0.3 


20.2±0.7 


3.6±0.5 


6.52 


16 


0.89 


1.09 


[W99]2 


I/U 


16.2±0.2 


6.4±1.0 


6.6±0.3 


14.2±0.4 


4.5±0.5 


6.31 


12 


1.54 


2.00 


[W99]l 


u 


17.5±1 


8.6±1.0 


8.1 ±0.5 


9.1±0.6 


3.6±0.5 


5.81 


19 


1.76 


2.22 


[W99]16 


u 


19±4 


9.9±1.0 


10±2 


15.8±1 


6.0±0.5 


6.51 


15 


0.50 


1.20 



ments is the reliability of the final results. For the optical 
echelle data, the instrumental resolution, (7instrument=-^-'2 
km s^^ ensures that the cluster line profiles are always 
well-resolved. The situation for the near-IR data is less 
clear. For the clusters where the velocity dispersions were 
estimated from the near-IR CO band-head, the measured 
values are above or approximately at the instrumental res- 
olution {crinstrument = ^4,.2 km s~^) and are thus likely to 
be secure. However, we were worried about possible sys- 
tematics effects that might influence the final estimates. 
To make an estimate of these possible problems, we con- 
structed "test data" from the stellar templates that were 
artificially broadened and had noise added to mimic the 
real cluster data from both the ISAAC and UVES obser- 
vations. This allowed us to investigate possible systematic 
problems in the way each of the data sets was analysed and 
allowed us to discover the most robust ways of analyzing 
the data. In addition, for one of the clusters, [W99]2, we 
obtained both UVES and ISAAC spectra of this cluster to 
observe any systematic differences between the two data 
sets that might lead us to conclude that there are system- 
atic differences in the properties among the clusters. 

From experiments with the test data, we found no sys- 
tematic problems in the estimated dispersions as long as 
the template was a good match (as described earlier) and 
for velocity dispersions that were above approximately 
50% of the instrumental resolution. We conducted these 
simulations for the range of S/N spanned by the data for 
the clusters (Figure ^). As a by product of these simu- 
lations, we were able to estimate the uncertainties in our 
estimates by conducting many (several tens) trials for each 
input broadening width. From the distribution of result- 
ing dispersion estimates, we were able to determine the 
average output velocity dispersion to look for systematic 



offsets and the variance to provide a likelihood estimate 
of the uncertainty in any one measurement. 

We applied the three criteria described in § [3.3.1 to 
our cluster data by creating a stellar template grid with 
CO(2.28/im) equivalent widths between 12. 2A and 18. 7A 
(we used the nine template stars listed in Table 2 and 
created intermediate spectra by averaging of two spec- 
tra), and added continuum levels of -30% to ±30%. One 
of these maps is displayed in Figure |[ The greyscale image 
of velocity dispersions has a contour overlay of the corre- 
sponding I and the minimum positions are to first order 
located at regions of comparable CO equivalent widths. 
We selected for the average a value only those template 
spectra which had a weak dependence of a on selected 
wavelength range, which is displayed on the second map 
in Figure ^, and additionally verified by visual inspection 
that the shape of the fit around the band-head is well- 
matched. 

The best- fitting template for cluster [W99]15, for ex- 
ample, is a star with a CO equivalent width of 17. 6A, with 
an additional 6% continuum contribution. This matches 
very closely the predictions of stellar contributions from 
evolutionary synthesis models: For example, according to 
STARS ([Sternberg 1998[ ), the K-band fight in a cluster of 
8.5 Myrs is composed of 95% M1/M2 type supergiants, 
~ 3% hot main sequence stars and wl.5% A-G type su- 
pergiants. 

The ISAAC and UVES observations gave an estimate 
of cr=14.0±0.8 and 14.3±0.5 km s"! for [W99]2, respec- 
tively, where the uncertainties in each measurement are 
the combination of the uncertainties in a single measure- 
ment as determined from the test data, and the deviations 
between different templates and fit wavelengths ranges. 
The difference in the estimated velocity dispersions is only 
0.3 km s~^. This difference is well within the random error 
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12.16 13.41 13.67 15.62 15.81 16.70 17.58 17.64 18.66 12.16 13.41 13.67 16.62 15.81 1670 17.68 17.64 18.65 

Wco(2.28nm) [A] Wco(2.28|.im) [A] 

Fig. 5. Visualization of the velocity dispersion determination for cluster W99-15. 

Left: Greyscale shows the velocity dispersion (black wlO kms~^ to white «27 kms^"'^), contours show the correspond- 
ing (from 1.28 to 2.0). Arranged along the x-axis are different input supergiant template stars, arranged according 
to increasing CO equivalent width. The input spectra between labelled columns were created by averaging two neigh- 
bouring spectra. The y-axis arranges the template spectra by added or subtracted continuum contribution, ranging 
from -30% to -1-30%. The line was added to guide the eye along the regions of lowest x^, which depends mainly on the 
CO equivalent width of the template spectrum matching that of the cluster. Some of the cr-values for the lowest ^ 
in a column are labelled, and the three best fitting a values are boxed, with the corresponding ^ value listed below 
the box. Right: The best fitting cr value is chosen from the three a values with the lowest (they are marked with 
boxes in the left figure) by requiring that a depend only weakly on the selected wavelength range. The corresponding 
diagnostic map is shown, which displays the standard deviation in the velocity dispersion (grey scale, ranging between 
0.05 (black) and 1.63 kms^^ (white)). Out of the pre-selected three best fits, the third spectrum, located between the 
columns labelled "17.58" and "17.64" shows the lowest variations in velocity dispersion for the five wavelength regions. 
The numbers in the boxes give the standard deviation for the velocity dispersions for the three best cases. Thus, the 
best matching template (marked by an ellipse) has a Wco of 17.6 A and a continuum contribution of 6%, resulting 
in a velocity dispersion of 20.2 kms^^, with a standard deviation of 0.37 kms^^. 



we have estimated using the test data and demonstrates 
that there are unlikely to be any systematic differences in 
using the UVES and ISAAC data (or the Tband compared 
to K-band or the CaT plus weaker lines compared to the 
CO band-heads) to estimate the velocity dispersions. 

We found that the best fitting models for the light pro- 
files of the clusters were King models with tidal to core 
radii ratios of 15 (c=1.176) or 30 (c=1.477; both concen- 
trations typically provided fits that were of indistinguish- 
ably good statistical significance but clearly better than 
the other fitted models, i.e., Gaussian or Moffat profiles). 
The effective radii of the observed clusters were approxi- 
mately 4 pc, similar to what Whitmore et al. (1999) esti- 
mated as the average radius of the population of clusters 
in the Antennae. So while the selected clusters are among 
the most luminous in the Antennae, they are certainly not 
among the largest. 



J^.l. Estimating the Metallicities and Ages of the Clusters 

Given the wide wavelength coverage of the UVES data, 
some interesting features are contained in the spectra. 
One such feature is the Mgl line at 8806. 8A. This fea- 



ture has been shown by Diaz et al. 1989 to be sensitive 



to the metallicity. The strength of this feature in two 
clusters is between 0.6 and 0.8A, which according to the 
empirical calibration of the equivalent width of the Mgl 
feature suggests approximately solar metallicity (there is 
significant scatter in the relationship between metallicity 
and line strength). However, the equivalent width of the 
Mgl absorption feature in [W99]2 is f«1.2A, which indi- 
cates super-solar metallicity. Unfortunately, for the clus- 
ters with only ISAAC near-infrared spectroscopy, it is not 
possible to estimate the metallicities of the clusters from 
these data alone. 
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Fig. 4. Artificial data with a S/N comparable to our data 
(S/N = 20) and with a S/N = 100 were created from steUar 
template spectra and the input velocity dispersion was 
redetermined. The average value of 30 test runs with the 
variance is displayed for both instruments (Top: ISAAC, 
Bottom: UVES). 

To determine the ages of the clusters, we followed the 
prescription outlined in Mengel et al. (2001), with the ad- 
ditional information from the CaT. This method relies on 
the equivalent widths of the CO-band head at 2.29/im, 
CaT lines, and the Br7 emission line (where available for 
each cluster). The Br7 emission strength estimates come 
from Mengel et al. (2001) or Mengel et al. (2001, in prepa- 
ration) with the later estimates resulting from narrow 
band Br7 imaging of the Antennae. We used the diag- 
nostic plots from Leitherer et al. (1999) of the equivalent 
widths of these features versus age assuming an instanta- 
neous burst, solar metallicity (but twice solar for cluster 
[W99]2), and a Salpeter IMF between 1 and 100 Mq. 
The strength of using these features is that the equivalent 
width versus age relations have different dependencies as 
a function of age. Br7 is strong in the early evolution of 
the cluster and weakens on the main sequence evolution- 



ary time scale of the most massive stars with significant 
ionizing photon output (ages less than about 6-8 million 
years). The CaT lines reach significant equivalent widths 
at about 6 million years (as the I-band light becomes dom- 
inated by red supergiant stars) while the CO band-head 
at 2.29/im becomes significant at about 8 million years (as 
the K-band light becomes dominated by red supergiants). 
The combination of these diagnostics have a useful and 
relatively unambiguous age sensitivity to instantaneous 
burst populations over the age range from to 20 mil- 
lion years, beyond which, the age determination ceases to 
be unambiguous. 

From these line strengths we estimated the ages given 
in Table 3, ranging from about 6 to 10 million years, which 
reflects the selection criterion of strong CO absorption, re- 
quired for the determination of the velocity dispersion (as 
discussed in § Moreover, we note that the metallic- 
ity estimates can play a crucial role in determining the 
ages. As we discussed for [W99]2, we have evidence that 
it has a super-solar metallicity. In Mengel et al. (2001), 
we found that [W99]2 exhibits strong CO band-heads as 
well as relatively strong Br7 emission and discussed the 
age determination of this cluster under the assumption 
of solar metallicity. Mengel et al. (2001) found that one 
age was difficult to reconcile with the measurements and 
speculated that the simultaneous Br7 and strong CO band 
head could be reconciled with a single age if the low spa- 
tial resolution of our data allow for different physical re- 
gions of different ages to contribute to the final spectrum, 
the cluster had an extended duration of star-formation, 
or that shocks were present. Our finding that this cluster 
may have super-solar metallicity releaves this dilemma in 
that it is now possible to reconcile all of the measurements 
in that strong CO absorption now appears at an earlier 
age (at an age where Br7 emission is also present). In 
this case, we now find the age of 6.6 Myrs instead of 8.1 
Myrs under the assumption of solar metallicity. Further 
we note, in our current sample, [W99]2 is that only clus- 
ter that shows significant Br7 emission and CO 2.29/xm 
lines simultaneously. 

5. Discussion and Implications 

5.1. Are the Clusters Bound and Have They Undergone 
Mass Segregation? 

The underlying assumption of any mass estimate using 
the velocity dispersion will be that the cluster is self- 
gravitating (i.e., bound). From the velocity dispersions 
and half-light radii we estimate that the crossing times 
{tcross ~ Thp/a) range from about 1 to a few xlO^ yrs. 
We determined the ages of these clusters to be around 8 
Myrs and thus these clusters have already survived for 20- 
50 crossing times. It is clear that at least initially, these 
clusters must have been gravitationally bound, and that 
the structure of these clusters is likely to be changing 
slowly over many crossing times. This assumption is also 
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supported by N-body simulations performed by Portegies 
Zwart et al. (1999) on the evolution of a cluster with a 
similar concentration like the ones in our study. Core ra- 
dius and density stay almost constant between 4 and 10 
Myrs, the endpoint of their study. Only one IMF was used 
(Scalo), but we assume that these findings should be valid 
unless the IMF slope is very shallow (below a « 1.5). 

The stars are expected to be segregated by mass on 
the "two-body relaxation time-scale" . Given that we have 
an estimate of the velocity dispersion and mass densities, 
it is possible to estimate the two-body relaxation time- 
scale. Using the arguments in Spitzer & Hart (1971) we 
estimate that the half-light relaxation time scale is several 
100 Myrs. Since the cluster ages are much less than this, 
dynamical mass segregation seems unlikely. However, sev- 
eral very young star clusters, like R136 in 30 Doradus (age 
around 4 Myrs, see e.g. Brandl et al., 1996) show evidence 
for mass segregation, which suggests that either the more 
massive stars formed preferentially near the center, or that 
the segregation occurs on a faster time-scale: N-body cal- 
culations by Portegies Zwart (2000) achieved some mass 
segregation within a few crossing times. If mass segrega- 
tion was present in the Antennae star clusters, both the 
half mass radius and the velocity dispersion, and therefore 
the mass, would be underestimated by our calculations. 

5. 2. The Estimates of the Masses 

Since we have argued that these clusters must have been 
gravitationally bound initially (although subsequent dy- 
namical evolution and mass loss through massive star 
winds and supernova explosions may unbind the clusters), 
we can now use the Virial theorem, the light profiles (and 
under the assumption that the light distribution follows 
the mass distribution, i.e., the mass-to-light ratio is not 
a function of radius), and the velocity dispersion to es- 
timate the masses of the clusters. The masses were de- 
rived assuming 1) that the best fitting profile was a King 
model (the light profiles were well fit with King models 
with concentrations of 1.176 or 1.477) 2) that the mass 
distribution follows the light profile, and 3) that the mea- 
sured velocity dispersions were in fact not the central ve- 
locity dispersions but were projected velocity dispersion 
weighted by the projected light profile. Again, the cor- 
rections between the half-light radius and the core radius 
(defined in the King model) and central velocity disper- 
sion and the measured velocity dispersion were given by 
those from the King models themselves (see, for exam- 
ple, McLaughlin 2000 for these relationships). The masses 
we derive are listed in Table 3 and range from about 6.5 
xlO^ Mq to 4.7x10*^ Mo). The uncertainties in the mass 
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estimates given in Table i.4 only represent the statisti- 



cal uncertainties estimated for the velocity dispersion and 
half-light radius and do not represent any additional un- 
certainties due systematic errors (as might be caused by 
a radial dependence on the mass-to-light ratio, system- 
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Fig. 6. The V-band light-to-mass ratios for the young 
compact clusters in the Antennae compared with models 
of Lcithcrer et al. (1999) for an instantaneous burst of to- 
tal mass 10^ Mq with a power-law slope, solar mctallicity , 
and a stellar mass range from 0.1 to 100 Mq. The power- 
law slopes are 2.35 (solid line), 2.00 (dashed line), and 
3.00 (dotted line). The points represent the data for the 
antennae clusters with V-band photometry from HST and 
with the extinction estimated as stated in the text. Un- 
fortunately, cluster [WS95]355 is too faint to be detected 
in the HST WFPC-2 images. The error bars represent a 
combination of the error in the mass and the uncertainty 
in the photometry and extinction estimate was assumed 
to be 10%. The light-to-mass ratios shown take into ac- 
count the mass lost through stellar winds and supernova 
over time. 



atically over-estimating the size of the cluster, erroneous 
functional form of the light profile, etc). 

We caution that while we believe the fits to the light 
profiles are robust, given that the HST data used to esti- 
mate the size scales has a resolution that is close to the 
measured half-light radii, it is possible that we have sys- 
tematically over-estimated the cluster sizes. This could 
lead to a systematic over-estimate of the total dynami- 
cal masses by the same proportion. The relation between 
cluster concentration and mass is non-linear. For our typi- 
cally best fitting models (c=1.176 or c=1.477), both lower 
or higher values would lead to a decrease in mass (by 2% 
if c=0.477 instead of 1.176, and by 25% if c=2 instead of 
1.477). This means that unless the clusters are much more 
concentrated than estimated, the masses can be overesti- 
mated only slightly. 



12 



Mengel et al.: Masses of Clusters in NGC 4038/4039 



5.3. Mass-to-Light Ratios of the Clusters and 

Implications for the IMF 

One of the critical parameters to understanding star- 
formation in and the dynamical evolution of these com- 
pact, young star-clusters is their initial mass function. 
If we find that these clusters are best described as hav- 
ing an IMF with significant numbers of low mass stars 
(say Salpeter slope down to 0.1 Mq), then the dynami- 
cal evolution will be driven through the ejection of these 
numerous low mass stars by binaries and massive stars in 
the cluster. Moreover, the fraction of the total mass lost 
through stellar winds and supernova explosion will be rel- 
atively small (10s of percent; Leitherer et al. 1999). It is 
then likely, since the amount of overall ejected mass is rel- 
atively small compared to the total mass, for the clusters 
to remain bound and long-lived. If however, there is a lack 
or relative deficit of low mass stars (either through an IMF 
that is truncated at the low mass limit or that has a rela- 
tively flat slope), then the total mass loss from the cluster 
due to stellar ejection (since higher mass stars would then 
be more likely to get ejected) and the fraction of the total 
mass lost through stellar winds and supernova would also 
be proportionally higher both of which would tend to in- 
crease the likelihood that cluster becomes unbound (e.g., 
Chernoff & Weinberg 1990; Takahashi & Portegies Zwart 
2000). 

In the following we use the derived masses and ages, 
in combination with the photometry and extinction esti- 
mates, to constrain the IMFs of these clusters through the 
use of population synthesis models (see for example, Stern- 
berg 1998). The approach uses a comparison of the cluster 
with population synthesis models to predict the mass-to- 
light ratios as a function of the cluster age. In our case, 
we can perform this comparison in both the V-band and 
the K-band, which is useful because these bands differ in 
their M/L-evolution with age and in their extinction sen- 
sitivities. For example, the K-band is less sensitive to the 
extinction than the V-band (by about a factor of 10) and 
also shows a strong dependence on age through the sig- 
niflcant effect on the K-band luminosity when supergiants 
become prevalent in the cluster stellar population. There 
are also practical differences in the two bands. The K-band 
data are ground-based and hence are sensitive to our abil- 
ity to remove the background from the photometry of the 
clusters and also might include contributions to their total 
magnitudes from other nearby clusters. The V-band data 
are HST WFPC-2 images and hence are less sensitive to 
contamination and background subtraction uncertainties. 
Although this last difference is partially mitigated by the 
fact that our K-band image used for this analysis has a 
seeing of «G'.'4 to G'.'5. 

Keeping these differences and limitations in mind, we 
attempt to constrain the functional form of the IMF in 

Figures 6 and 7. In the plot of the log Lk /M versus age, 
the most striking feature is that the clusters form two 
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Fig. 7. The K-band light-to-mass ratios for the young 
compact clusters in the Antennae compared with mod- 
els of Leitherer et al. (1999) for an instantaneous burst of 
total mass 10^ Mq with a power-law slope, solar metal- 
licity (we did not add the track for 2xZq, which would 
be appropriate for cluster [W99]2because the difference is 
below 0.1 dex, well within our uncertainties), and a mass 
range from 0.1 to 100 Mq. The power-law slopes are 2.35 
(sohd line), 2.00 (dashed line), and 3.00 (dotted hue). The 
points represent the data for the Antennae clusters with 
ground-based K-band photometry (Mengel et al. 2001, in 
preparation) and with the extinction estimated as stated 
in the text. The error bars represent a combination of the 
error in the mass and the uncertainty in the photome- 
try and extinction estimate was assumed to be 10%. The 
model light-to-mass ratios shown take into account the 
mass lost through stellar winds and supernovae over time. 



groups, several clusters have Lk /M at their estimated ages 
that seem best described with an IMF having a relatively 
flat slope, or perhaps slightly steeper slopes, but truncated 
low mass limits (higher than the 0.1 Mq in the model). 
For example, the difference in mass between a Salpeter 
IMF slope with an initial mass range from 0.1 - 100 Mq 
compared to a Salpeter IMF slope with an initial mass 
range from 1.0 - 100 Mq is a factor of 2.6. Such a differ- 
ence would shift the Salpeter IMF slope line in Figures 6 
and 7 upwards by 0.4 dex implying that a Salpeter IMF 
slope with an initial mass range from 1.0 - 100 Mq is con- 
sistent with the results for clusters [W99]l and [W99]2. 
The results for clusters [W99]15 and [WS95]355 are more 
consistent with a steeper IMF slope, approximately 2.5 
with an initial stellar mass range of 0.1 - 100 Mq or a 
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steeper slope with a truncated lower mass. In the light- 
to-mass ratios in the V-band, this difference is less clear. 
The light-to- mass ratios for [W99]l and [W99]15 give con- 
sistent results in both the K- and V-bands. 

Similar results for young, compact clusters in nearby 
galaxies have also been presented. Sternberg (1998) us- 
ing data from the literature for two clusters, one in NGC 
1705 (NGC1705-1) and another in NGC 1569 (NGC1569- 
A), found that the light-to- mass ratio of these two clusters 
indicated rather shallow IMF slopes (« 2.5 for NGC1569- 
A and < 2 for NGC1705-1 if their IMFs extend down to 
0.1 M0). Smith & Gallagher (2001) found in their analy- 
sis of M82-F that it either requires an extremely flat IMF 
slope (a << 2.0) or that it has a steeper IMF slope (a — 
2.3) but then has a lower mass cut-off of about 2-3 Mq. 
However, some caution is necessary in the interpretation 
of M82-F. The photometry of this cluster is rather un- 
certain given that it was heavily saturated in the HST 
images used by Smith & Gallagher (2001) to estimate its 
magnitude and the extinction (Es-y =0.9±0.1). Smith 
& Gallagher (2001) also re-analysed some of the half-light 
radius and magnitude measurements of NGC 1705-1 and 
NGC 1569- A and although the detailed estimates of the 
light-to-mass ratios for NGC1705-1 differ, they come to 
the same conclusions as Sternberg (1998). We show a 
graphical comparison of our results with those obtained 
for M82-F, NGC1705-1, and NGC1569-A from Smith & 
Gallagher (2001) in Fig. 8. 

If the segregation of the cluster points in Figure 7 is 
a real effect, it means that the IMF varies significantly 
between clusters. We do not yet have the possibility to 
disentangle lower mass cutoff and IMF slope, but the two 
different groups of clusters differ in their ratio of low- to 
high mass stars by approximately a factor 5 ("low" and 
"high" relating to stars below and above 3 Mq, respec- 
tively) . 

There has been substantial debate in the literature as 
to whether or not different environments lead to different 
forms of the IMF. While there does not exist a thorough 
understanding of the process of star-formation, possible 
mechanisms that have been investigated for determining 
the characteristics of the IMF are variations in the local 
thermal Jeans mass (e.g., Larson 1985; Klessen & Burkert 
2001), turbulence (e.g., Elmegreen 1997), gravitationally- 
driven fragmentation (e.g., Larson 1985), cloud-cloud and 
star-star coalescence (e.g., Nakano 1966; Bonnell, Bate, & 
Zinnecker 1998) and several other processes (see for exam- 
ple, Larson 1999| , Elmegreen 1999 and references therein). 
While it would be highly speculative to suggest that the 
range of IMF slopes and/or low mass cut-offs we observe 
supports or refutes any of these hypotheses for what phys- 
ical processes are responsible for determining the charac- 
teristics of the IMF, it is useful to discuss what constraints 
these results might offer in understanding the origin of the 
IMF. 
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Fig. 8. The points are the same as in Fig. ^but now with 
additional results from Smith & Gallagher (2001) for the 
clusters M82-F, NGC1705-1, and NGC1569-A. The results 
obtained for the clusters NGC1705-1 and NGC1569-A are 
similar to those in Sternberg (1998) except for NGC1705- 
1 which in Smith & Gallagher (2001) has a Ly/M lower 
by a factor of approximately two mainly due to their re- 
evaluation of the cluster half-light radius (revising the size 
upward by about a factor of 2). 



There are two attractive hypotheses for explaining the 
observed segregation. One is that perhaps the difference 
in the relative number of low mass stars is related to 
the environment in which the cluster is born. Clusters 
[W99]15 and [WS95]355 suffer greater extinction on av- 
erage than most clusters in the Antennae and both lie 
in the over-lap region of the two merging disks. Clusters 
[W99]l and [W99]2 lie in the "western loop" of the Anten- 
nae's northern most galaxy, NGC4038 and have relatively 
low extinctions (AyftiO.3 magnitudes). In spite of the clus- 
ters having similar ages, the environments of [W99]15 and 
[WS95]355 appear to be much more rich in dust and gas 
(Whitmore & Schweizer 1995; Wilson et al. 2000; Mengel 
et al. 2001). However, cluster [W99]16, while physically ly- 
ing near the clusters [W99]15 and [WS95]355, appears to 
suffer much less extinction. Presumably, since the region 
where [W99]15, [W99]16, and [WS95]355 were formed is 
also the region where the two disks are interacting, the 
gas pressure (both the thermal and turbulent) would be 
high due to dynamical interactions of the two interstellar 
media. The local thermal Jeans mass, which is considered 
to play a role in setting the mass scale for low mass stars 
(e.g., Larson 1999), is proportional to x^P^^/^, where T 
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and P are the total pressure and temperature (e.g., Spitzer 
1978). Clouds with low temperature and/or high pres- 
sures would form clusters with IMFs weighted towards low 
mass stars, which would be applicable for the more extin- 
guished clusters. Conversely, clouds that are collapsing in 
relatively low extinction, gas poor, more quiescent regions 
might be influenced by the ionizing radiation of nearby 
young clusters raising its gas temperature. Such a situa- 
tion would favor the formation of clusters with IMFs that 
have relatively large fractions of high mass stars. Such an 
idea might also explain the mass functions of NGC1569-A 
and NGC1705-1 whose relatively low metallicities could 
favor high gas temperatures (by reduced cooling and high 
background radiation fields). Thus studies of young com- 
pact clusters in the Antennae and other galaxies appear 
to agree with the expectation from the low mass scaling 
set by the Jeans mass. 

The second plausible hypothesis is that in Fig. 7 and 
8 the clusters with relatively large ratios of light-to-mass 
are also among the lower mass clusters. N1569-A, N1705- 
1, M82-F, [W99]l, and [W99]2 all have masses less than 
about 2 X 10^ Mq. Clusters [W99]15, [W99]16, [WS95]355 
all have masses greater than 3 x 10^ Mq. While the total 
mass of a cluster influencing the observed IMF has been 
discussed in the literature, it is most commonly discussed 
within the context of random sampling of stars in clusters 
with relatively few stars (e.g., Elmegreen 1999). The An- 
tennae clusters are certainly not in this limit since they 
have many x 10^ stars. However, it may be that there is a 
process internal to the cluster that inhibits preferentially 
the formation of low mass stars and that this may scale 
roughly with the mass of the cluster. One possibility is 
that the massive stars through radiation pressure might 
inhibit or shut off the formation of low mass stars. 

Obviously using compact, young star clusters in the 
Antennae and other galaxies to provide observational con- 
straints on the processes that lead to the characteristics of 
the IMF will require much more observational and theoret- 
ical work. However, even the relatively sparse and modest 
results provided here and in the literature are tantalizing 
and show the promise of this type of observation in ad- 
vancing our understanding of the process of cluster and 
star- formation. 

5.4- The Evolution of the Most Massive Clusters: Will 
They Dissolve? 

Are these young clusters progenitors of the old globu- 
lar cluster systems that we observe in normal galaxies? 
If they are, how would the population of young clusters 
need to evolve to have a mass distribution that resembles 
that of the old globular cluster system? These fundamen- 
tal questions can only be addressed when we have char- 
acterized the masses, initial mass functions, and strength 
and change in the gravitational potential in which these 
clusters form and evolve. 



Various models have been constructed to fol- 
low the evolution of clusters with structure simi- 



lar to globular clusters (see Spitzer 1987 and refer- 
ences therein). Perhaps the most relevant for this 
study are those that employ the Fokker-Planck ap- 
proximation in solving the dynamics of the constituent 
cluster stars (|Cohn 1979|, phcrnoff fc Weinberg 1990 



Takahashi fc Portegies Zwart 200C| )! The dynamics of 
Fokker-Planck models are relatively simple to calculate 
compared to full N-body models and thus other as- 
trophysical effects that might be important to the dy- 
namical evolution of the cluster can be included (such 
as including a range of stellar masses distributed with 
a realistic IMF, stellar mass loss, stellar evolution, 
etc). While these models have their limitations, the re- 
sults of these studies provides a useful guide on how 
these young, compact clusters of the Antennae might 
evolve given our determined constraints on the shape 
of their IMFs and concentrations. The Fokker-Planck 



models (e.g., 


Cohn 1979, 


Chernoff & 


Weinberg 


1990 


Takahashi & Portegies Zwart 


200C|) sug 


gest the 


most 



likely survivors are those clusters that are either more 
highly concentrated (King model concentrations of greater 
than 1.0-1.5), and/or those with steep IMFs. This is shown 
graphically in Fig. 9 where we have reproduced approx- 



imately a figure from Takahashi & Portegies Zwart 2000 
(their Figure 8) that shows which range of cluster param- 
eters the clusters are long-lived or disrupt. For the range of 
plausible IMF slopes that we have determined (assuming, 
as in the models, that the IMF extends down to sub-solar 
masses), the models suggest that clusters that are about 
as concentrated as we have observed in the Antennae will 
survive for at least a few Gyrs. Even though we have no 
independent constraint on the lower mass limit of the stel- 
lar mass function, it is interesting that models with steep 
IMF slopes are also robust against disruption. From these 
models alone it is difficult to tell if the clusters are able to 
withstand the disruptive tendencies of the large mass loss 
they suffer if the lower mass limit (0.1 Mq in Figs. 6 and 
7) is raised and the IMF slope steepened, such that the 
observed light-to-mass ratio constraints are met. The lim- 
ited number of models so far generated seems to indicate 
that clusters with very steep slopes {a > 3) of their IMF 
are robust. More modeling is obviously needed to answer 
this question quantitatively and definitively. 

Even if the clusters survive, it is obviously impor- 
tant to have some indication of how much mass loss 
they may undergo as they age. For some of the clus- 
ter parameters, these models suggest that these clus- 
ters may lose significant amounts of mass (>50%, some- 
times much greater; Takahashi & Portegies Zwart 2000). 
Clusters with steep IMF slopes {a > 2.5) which 
are compact (c>1.5) apparently suffer large mass 
loss (may lose more than 90% of their total ini- 
tial mass; [Takahashi fc Portegies Zwart 200C but also 



[Chernoff fc Weinberg 199C whose models suggest that 
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Fig. 9. This is a modified reproduction of a figure from 
Takahashi fc Portegies Zwart (2000; their Figure 8) show- 
ing the ages at the endpoint of evolution for modeled clus- 
ters. The clusters were characterized by scaled central po- 
tentials, Wo (which is related to the concentration, c, as 
shown in the ordinate on the right side) and IMF slopes 
(q). The numbers in the grid give the ages at the endpoint 
of the evolution (disruption or core collapse) as log(age 
[Grys]). Clusters with values below 0.3 (corresponding to 
2 Gyrs) were disrupted by dynamical processes over time, 
the others experienced core collapse. The grey box shows 
the region where the clusters in the Antennae are observed 
to lie. This suggests that the clusters are relatively long 
lived. 



perhaps some of these high mass loss clusters are dis- 
rupted). Therefore, it is easily possible for these cluster 
to survive, but over their lifetimes, they will have much 
lower masses than when they are as young as the clus- 
ters we have observed here. It is through such strong mass 
loss that the young compact clusters in the Antennae may 
evolve into something similar to the population of clusters 
observed in galaxies generally (e.g., Pritze-v. Alvensleben 
1998; 1999). 

There are limitations in the models that might in fact 
lead to higher mass loss rates and possibly easier disrup- 
tion of clusters for a given set of initial conditions than 
predicted in the models. One is that the tidal field is much 
more complex and stronger than what is assumed in the 
models which is static or a simple orbital potential in a 
disk. In a strong and variable tidal field, the mass loss 
rates would undoubtedly go up (e.g., Gnedin & Ostriker 
1997). However, very little is known about the tidal field 



over time on the scales of the clusters and such an analy- 
sis would require a greater understanding of the mass dis- 
tribution in galaxies than we currently have. Therefore, 
literal interpretations of the evolution of clusters as indi- 
cated by the models should be viewed with some caution. 

In summary, given what we have learned in this anal- 
ysis, it is likely that the most massive clusters arc able 
to survive in the galaxian environment in which they are 
born, and do not abruptly disappear. Since we cannot con- 
strain the lower mass limit of the stellar mass function it 
is possible that the clusters are easier to disrupt than we 
have suggested here. The amount of mass- loss they will 
undergo is uncertain but is likely to be a large percent- 
age of their initial masses thereby providing a mechanism 
were the massive clusters studied here can evolve to look 
more like a typical cluster in a galaxy like the Milky Way. 

6. Summary and Conclusions 

In order to estimate the masses of young compact clusters 
in the interacting pair of galaxies NGC 4038/4039 ("the 
Antennae" ) , we have measured the velocity dispersions of 
6 clusters using high resolution optical and near-IR spec- 
troscopy conducted at the ESO-VLT. The velocity disper- 
sions were estimated using the stellar absorption features 
of CO at 2.29 /um and metal absorption lines at around 
8500 A including lines of the Calcium Triplet. To estimate 
the masses, we measured the size scales and light profiles 
from archival HST WFPC-2 (see Whitmore et al. 1999). 
The ages of individual clusters were estimated using the 
equivalent widths of the CO band- head at 2.29 /im, the 
CaT, and the Br7 emission line (estimated from data in 
Mengel et al. 2001 and Mengel et al. 2001, in preparation). 
Our principle results from this analysis are: 

— The clusters have ages of about 8 Myrs. This is a selec- 
tion effect since wc selected the clusters based on other 
observations indicating that they had ages of about 10 
Myrs to ensure they had strong COA2.29 jim and CaT 
absorption lines for determining the velocity disper- 
sions. 

— The best fitting light profiles are consistent with King 

models of concentration parameters, c, of either 1.176 
or 1.477 (tidal to core radius ratios of 15 or 30 re- 
spectively). The half-light radii from these fits are all 
around 4 pc which is typical of the young compact 
clusters in the Antennae (Whitmore et al. 1999). 

— The measured velocity dispersions were all between 9 
and 21 km s""'^. Combining the velocity dispersions, 
half-light radii, and parameters of a King model under 
the assumption that the clusters are in approximate 
Virial equilibrium (an assumption which is supported 
by our estimate that the clusters have already survived 
for 10s of crossing times) yields mass estimates that 
range from 6.5 xlO^ to 4.7 xIO^Mq. These masses 
are large compared to typical masses of a globular clus- 
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ter (~ 1 X IO^Mq), but this not surprising, since we 
selected the brightest clusters for our first analysis. 

— Comparing the cluster light-to-mass ratios with stellar 
synthesis models (Leitherer et al. 1999) suggests that 
these clusters may exhibit a range of initial mass func- 
tions, with clusters in the "over-lap" region showing 
evidence for a steeper IMF slope than those clusters 
in less extinguished regions. However, with only five 
clusters observed this result needs to be substantiated 
with observations of more clusters. 

— Results of Fokker-Planck simulations of compact clus- 
ters with concentrations and IMF parameters similar 
to those found for the young, compact clusters in the 
Antennae studied here, suggest that these clusters are 
likely to be long-lived (more than a few Gyrs) and may 
lose a substantial fraction of their total mass. 
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